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Abstract 

We consider a one dimensional evolution problem modelling the dynamics of an acous- 
tic field coupled with a set of mechanical oscillators. We analyse solutions of the system 
of ordinary and partial differential equations with time-dependent boundary conditions 
describing the evolution in the limit of a continuous distribution of oscillators. 
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1 Introduction 

We describe first the physical setup our system of equations is meant to model. We consider 
the following arrangement: an infinite pipe filled with a nonviscous, compressible fluid and 
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n mechanical oscillators realised with thin walls of mass Mj positioned in the pipe perpen- 
dicularly to its axis. A system of springs with elastic constants Kj confine the walls around 
their equilibrium positions located in the points of coordinates Sj, with respect to a system of 
coordinates whose x-axis coincides with the axis of the pipe. 

We suppose that there is no friction between the fluid and the pipe and we analyse only 
one dimensional flows. The acoustic field is then described by the pressure field p(x, t) and 
the velocity field v(x,t). The motion of the mechanical oscillators is described through the 
displacements yj(t) of the thin walls from their equilibrium positions Sj and the velocities 

z j(t) =%•(*)> j = 1, 

The field p(x, t) specifies the deviation of the pressure in the point x at time t with respect to 
an equilibrium pressure P . We analyse the dynamics of the whole system (acoustic field and 
mechanical oscillators) in the linearized acoustic and elastic regimes. Moreover we consider 
walls of zero thickness and we suppose that the contact between the fluid and the walls is 
maintained throughout the evolution. 

Under these hypotheses the dynamics is described by the following system of ordinary and 
partial differential equations with time dependent boundary conditions 



x e R\S (2) 



,f = -aV - * eR \s (i; 

dv 1 dp 
dt p Q dx 
dz- 

= -K m - S (p( a +) - p(sj)) j = 1, . . . , n (3) 

*i = = J = l,...,n (4) 

where po is a positive constant representing the density of the fluid, a is a positive constant 
indicating the velocity of sound in the fluid, S is the area of the transversal section of the pipe 
and S := (si, . . . , s n ) is the set of the equilibrium positions of the walls. 

Aim of this paper is to characterise the effective behaviour of the system in the limit of a 
continuous distribution of oscillators in a bounded region. 
Given the initial conditions 



p(x, 0); v(x, 0); %(0); Zj(0) j = 1, . . . , n (5) 
we consider the Cauchy problem ((!])- (jlj), Qj. 

As it was shown in [I] existence and unicity, together with detailed properties of the solutions, 
can be proved formulating the problem as a unitary evolution in the space of finite energy. In 
the following we fix the notation and we recall the main result obtained in [I]. 
We indicate with a capital Greek letter the generic vector (p, v,y,z) E L 2 (R) ©L 2 (1R) ©C n ©C n , 
where L 2 (M) is the space of square-integrable functions on the real line, 



V = 2/1^1 H h y n e n , z = z x e x H V z n e, 

and e 1; . . . ,e n is the canonical orthonormal base in C n . 
The Hilbert space of finite energy is defined by 

U := L 2 {R) © L 2 (R) © C n © C n 
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with the scalar product 



1 1 n 

((# i, * 2 )> := -3— (Pi,P2) + Po(ui, ^2) + ^ V [i^- j/yj/2j + Mj zy^y] , 
a po '■ * 



3=1 



where (• , •) indicates the standard scalar product in L 2 (R) and ~~ denotes complex conjugation. 
The square norm of a vector \1>, ||^|| 2 = ^)), defines the total energy of the system in the 
state * 



^toi = 77 II ^11 2 = Eac + E osc , 



with 



EaC ~ 2a 2 Po 



J -00 z J -c 



\v(x)\ dx 



E = V 

osc / 
3=1 



E'ac is the energy stored in the acoustic field while i£ osc is the energy of the mechanical oscillators. 
Let H 1 (R) be the homo geneous Sobolev space of locally square-integrable functions with square- 
integrable (distributional) derivative and if-^R) the usual Sobolev space if 1 (R) := H 1 (W) fl 
L 2 (R). Similarly we define, for any open set OCR, the Sobolev spaces H\0) := {/ G L 2 (0) : 
£ G L 2 (0)} and H 2 (0) := {/ G if^O) : G L 2 (C)}. The following theorem holds 



Theorem 1. The linear operator 



A : D(A) CH^H 



D(A) := {(p,v,y,z)\pe L 2 (R) n H\R\S), v G H 1 ^), 
y G C n , z G C\ u(sj) = Zj, j = l,...,n} 

dv 1 cip re 



is real and skew- adjoint. Here Oj G C is defined by 



, ' m ± up req \ ~\ / if,- 5 



and p reg G f/" 1 (R) zs £/ie regular part ofp, i.e. 



1 n 

p res (x) := - - a 3 s S n - s i) 



i=i 



Remark 1. Notice that the left and right limits entering in the definitions of the Oj's exist and 
are finite by 



\p(xi) -p{x 2 )\ < (x 2 - x x 



.1/2 



dp 



, Sj < X\ < x 2 < s 



L 2 ( Sj ,s j+1 ) 



3+1 
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Theorem [T] implies that the Cauchy problem 



±9 t = A9 t (6) 

% =0 = V en (7) 

has a unique solution expressed via the strongly continuous unitary group of evolution generated 
by A, ty t = e tA ^Q. Moreover the unitary evolution preserves reality in such a way that one can 
consider the flow restricted to the real Hilbert space 



u r ■■= mm ©Li 



in 



1 



where L 2 r {R) is the space of real- valued functions in L 2 (R). The operator A is skew-adjoint so 
that H^tll = || ||, and the total energy remains constant. It is easy to check that the unique 
solution of (E])-(J7]) solves in H the Cauchy problem (JTJ)- (SD , ©• 
In fact from Theorem [1] the equation ([6]) is equivalent to the system 

dp 2 

dt a ^° dx ' 



dv 1 dp. 



reg 

> 



Zj , j 1, . . . , n . 



dt po dx 
dyj 
dt 

dzj ( Kj S 

Notice that if \l/o £ D(A) the solution ty t belongs to D{A) so that condition (j4j) is automatically 
satisfied at any time. 

In the following section we investigate the effective response to an incoming acoustical wave of 
a chain of mechanical oscillators in the limit of a continuous distribution of oscillator masses. 



2 Effective equation in the limit of a continuous distri- 
bution of oscillators 

Let us consider a sequence of Cauchy problems 

dt 1 1 

for increasing n G N, when the following assumptions are satisfied: 

A\) The set := (s^\ . . . , s^) of the equilibrium positions is contained in a bounded 
interval around the origin, Sj G [-L/2, L/2] Vn and \fj = 1 . . . n; 

Aq) The two measures fi^ := Y^j=i Mj^5 a (n) and [i^ := Y^j=i ^j"^ s C™)> where 5 s ( n ) is the 
Dirac mass supported at point are uniformly equivalent in the sense that < c\ < 
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sup [L 2 Kj U ^ I a 2 < C2 < oo, here c\ and C2 are two positive constants. The two mea- 
sures have total mass which is bounded uniformly in n : Y^j=i < c Y^j=i Kj < c - 
Furthermore and p^ weakly converge, when n tends to infinity, to measures sup- 
ported in [-L/2, L/2] which are absolutely continuous with respect to Lebesgue measure 
(the corresponding Radon-Nikodim derivatives will be denoted with SpM and Spx)] 

A3) Initial conditions will be chosen independent of n. In particular we consider the case 

p ( - n \x,0)=p (x); 0) =!*(*); ^ (n) (0) = 0; zf(0) = j = l,...,n, (8) 

where p and v are in H 2 (R) and supp(p ) n [- L/2, L/2] = 0, supp(w ) n [-L/2, L/2] = 0. 

Remark 2. The request of equivalence of the two measures p^ and p£ could be released in 
various ways but it is well motivated, from a physical point of view, preventing the occurrence 
of zero or infinite proper oscillation frequencies. 

Initial conditions more general than (jSJ) might of course be considered. Nevertheless OH]) seem to 
be the most natural initial conditions which do not depend on n. Notice, in particular, that at 
each step n, corresponding to the zero eigenvalue of A^ n \ there are static solutions of problem 
(l)-(4). As it was shown in [1] the static solutions are of the form fy st = (p Jt , 0, y , 0) with 
supp(p s t) C [—L/2, L/2]. Being initial conditions (jHJ) orthogonal to all static solutions we are 
guaranteed that a genuine scattering problem is analysed. 
Under the assumptions stated above we prove the following 

Theorem 2. Let the sequence of sets := (s[ , . . . , Sn) and the sets of positive real numbers 
M (n) := (M[ n) , M ( n n) ) and K {n) := (K[ n) , K { n n) ) be such that the assumptions A x ) and 
A2) above are satisfied. Let v^ n ' be the velocity field in the solution of problem (25)- (26) with 
initial conditions satisfying assumption A3). Then there exist a subsequence 

c {v {n) }™ C C 2 (R; L 2 (R)) n C\R; H\R)) n C(R; H 2 (R\S (n) )) 

and 

v e C^R; C(R)) n C(R; H^R)) 

such that 

1. t>( nfc ) and 9v< Q t k) converge to v and ^ uniformly over compact subsets ofR 2 ; 

converges to ^{t, •) in L 2 oc (R), uniformly in t over compact intervals; 
3. v is a weak solution of the hyperbolic equation 

, , Pm\ d 2 v 2 d 2 v p K , . 



with initial data 

dv 

v\ t =o = v ; - 



t=o 



1 dp Q 
p dx 



(10) 



Proof. By the basic properties of one-parameter groups of operators and their generators (see 
e.g. 0), the map t ^ e' i(n) ^ belongs to C^R; H)nC(R; D(A^)) if and only if ^ e £>(i (n) ). 
In this case 

dt 



and, since A^> is skew-adjoint, 



(11) 



Therefore, since any initial condition \l>o = (po,^o,0 ? 0) chosen according to assumption A3) 
is in D([A^] 2 ) for any n, the path t H> ^[ n) := e tA(n) ^/ is in C 2 (M;H) n C\R; D(A^)) n 
C(K;£([i (n) ] 2 )), 



dt 2 



(12) 



and 

||[i(n)]2^W|| = ||e^ (n) [iH] 2 ^ || = ||[i(")] 2 * ||. 
From the explicit characterization of domain and action of [A^] 2 



(13) 



D{[A in) ) 2 ) = \ [p (n \v {n \y^ n \z {n) ) E D(A {n) ) 
dp ™ 9 EH\R), ^(s^)= Po 



dx 



dx 3 



K (n) 



(n) y J 



y ' + 



dv {n) 
dx 

S 



e H l {R\S {n) ), 



» 



(n) C 3 ' ■ 



j = l,...,n 



[A^}\ P ^M n) J n \^) = ( a 2 ^, « 2 ^(^ 



y^( j .. < " > 



i\M : 



(n) y j 



reg 



(n)°i ' 



J VM, 



(n) 



M, 



(n) ^ 



where 



(n) dt) (n) (n)+ (n) _ 
C := — — is, ) : — (s) 



dx J 



dx 



of := V ^\sf^)-^\sf- 



dx 



dv {n \ , 
{x):= cJx- {x > 2 

reg j =1 



n 

1 An) I 



(14) 



We conclude that if ^J n) = (p^ (t) , (t) , (t) , z}^ (t)) is the solution of flT2}, then 



,(») G C 2 (K; L 2( R ^ n C 1( R . #1( R )) n ^(R; ff2(R\5(«))) , 2 g c 2 



and n )) is the strong solution of the following Cauchy problem with time- dependent 
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boundary conditions: 



dt 2 <).r \ Ox n (/ 

j2 (n) 7^(n) „2,. o 



(«) 



dt 



,(n) 



li=0 







1 <9po 
= 0. 

t=0 



(15) 
(16) 

(17) 
(18) 

(19) 



Furthermore, by energy conservation and by f ill I) and (|13[) . there exists a positive constant c 
depending only on the norm of the initial datum ty = (po, vo, Q, 0), such that the following 
bounds hold for any t and n: 



v (n) {t,-)\ 



H 1 (R) 



< C. 



Q v (n) 



Ot 

Q v (n) 

dx 



< c. 



H 1 (R) 



dt 2 



L 2 (R) 



reg 



<C. 



K 



(n) 

i _(«) 

(n) Z j 



a 2 poS A ( n ) 



i=i i=i 
By the continuous injection if 1 (IR) C&(R) one has 

,(») 



(t,a;)eR 2 

and (here Xi < x 2 and tj. < t 2 ) 

\v ( - n \t 1 ,x 1 )-v^(t 2 ,x 2 )\ 

Q v (n) 



sup |*/ nj (t,x)| < cisup||v (,l) (^-) 



lfr!(M) 



< c. 



(20) 
(21) 
(22) 



<\fx 2 — X\ SUp 
tgR 



<9x 



+ (*2— tl) SUp 

L 2 (R) (t,x)€B 



<y/x 2 ~ x[ SUp ||w (n) (^, -)||hi(R) + Cl ^ 2 ~ *0 SU P 



dt 

Q v {n) 



(t,x) 



dt 



(t, 



i3"!(R) 



Hence by (!20|) the sequence is equi-bounded and equi- continuous. Thus, by Ascoli- 

Arzela theorem, there exist a subsequence {i/ nfe )}^° and v E C(IR 2 ) such that ?/ nfe ) converges to 
v uniformly over compact subsets of M 2 . 

Let us now take a bounded rectangle R — I\X I 2 C K 2 . Since the injection H 1 (I 2 ) > C{I 2 ) 



is a compact operator and, by ( 120|) . 



9t 



is bounded in L°°(Ji; H l {I 2 )) and 



at 2 



is bounded in L°°(Ji; L 2 (J 2 )), by corollary 4 in j3], section 8, applied to the triple of Banach 
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spaces X = H 1 {I 2 ) C B = C{I 2 ) C Y = L 2 (I 2 ), there exists a sub-subsequence | 9v q^ j 
converging to vr in C(/i; C(J 2 )). 

Let H a (I 2 ), a G (0, 1), be defined, as usual, by 

H a {h):={feL\l 2 ) : / / ^ dxdy < oo 



h Jh F - 1/1 



H a (l2) is a Banach space with norm 




h 2 Ji 2 \x-V\^ M J 

Notice that sgn (■ — Sj) G H a (l2) for any a < 1/2 and for any bounded I 2 C R and its H a (I 2 ) 
norm is bounded uniformly in n. Hence, by (fill) . (12"T|) . (I2"2l and .A 2 ), { 9v( " k) X is bounded in 



I 9 * Ji 

L°°(Ii;H a (I 2 )), a < 1/2. Moreover, by flU}, {|^^}~ is bounded in L°°(h; L 2 (I 2 )). Thus, 

since the injection H a (I 2 ) L 2 (I 2 ), a > 0, is a compact operator, again by corollary 4 in [3], 
section 8, applied to the triple of Banach spaces X = H a (I 2 ) C B = Y = L 2 (I 2 ), there exists 

f ^ jj ^ CO 

a sub-subsequence < dv Qx > converging to v' R in C(I\\ L 2 (I 2 )). 

Considering then the sequence of rectangles R m = (— m, m) x (— m, m), by a standard 
diagonal argument, one obtains a further subsequence, which for notational convenience we 

continue to denote by {t> (n ' : )}^ , such that j 9v ^ - | converges to w G C(R 2 ) uniformly over 

compact sets, and such that j 9v ^ x k) j converges to i>' G C(R; L 2 oc (R)) uniformly over compact 
time intervals. Defining 

ft PX 

Vi(t, x) := v (0, x) + sgn (t) I v(s,x)ds, v 2 (t, x) := v(t, 0) + sgn (x) / v'(t,y)dy, 

Jo Jo 

one obtains in"*) — >■ t>i and t>(™ fe ) — >■ t> 2 , and so 

<9t> . <9t> 
wi=i; 2 = w, v = — , v = — . 

at ax 

Let us now multiply both sides of equation (fT51) by a function <? G C^R 2 ). Integrating by 
part we get 

d v {n k ) ^ s dg fJ _ ,\ u /dv^ nk \ n , _ /n x \ 2 f°°/dv^^ s dg 



x a (v),^(«,o;*-^- s -(o,-), 9 (o,-); = -a'jf ^-ar(v),^(v))«« 



Using assumption A2) and .A3) we conclude that 

1 / ^nAn k ) aZ j f^Og 



k^oo p S J dt y ' dt 1 

lim J- fVM^ r^rt -^(t 



and 



I />oo n fc 

lim — — — / V K<p k) zf k \t)g(t, sf k) ) dt 

= , lim -o / E " # fc) )) 9(t, *t k) ) dt 

1 r°° n 

+ , lim / E 4* fc) 4 nfc) ) * 

= J o °°(^v(t,.),g(t,.)\ dt. 

The limit function t> then satisfies the equation (notice that §|(0, •) and pm have disjoint 
supports) 

-f (( 1+ ^)l ft - ) -l ft - ) >*-(( 1+ ^)i (R - ) -> 



-a 2 / *( t ,.)\ A- r(EE v(t ,.), g{ t,.) 



o \dx 7 dx ' / Jo \ Po 
which is a weak form of equation fl9]) when the initial data v is specified. □ 
Remark 3. Posing 

1 . PM PK 

w:=l-\ , q:= — , 

Po Po 

let us define the weighted L 2 space 

L 2 {R; w) := {/ : R ->■ C measurable : /v 7 ^ G L 2 (R)} 
and the linear operator 

L : £>(L) C L 2 (R; w) -> L 2 (M; w) , L/ := — f-a 2 ^ + ?/ 



D(L) := If G L 2 (M; u>) : / G C^R) , ^ G AC(R) , L/ G L 2 (R; w] 



Here AC(R) denotes the set of absolutely continuous functions on the real line. Notice that L 
is in the limit point case at both +00 and —00 (use e.g. Theorem 6.3 in [4]) and so, by the 
theory of Sturm-Liouville operators (see e.g. [I]), L is self-adjoint (and positive). Hence, by 
the theory of abstract second-order equations (see e.g. [2], Chapter 2, Section 7), the Cauchy 
problem 

d 2 v 2 d 2 v 



dt 2 dx 2 

dv 

v\t=o = vi G D(L) 



= v 2 e D(VZ) 

t=0 



dt 

has an unique strong solution 

v G C 2 (R; L 2 (R; w)) D C^R; D(VZ)) n C(R; D(L)) . (23) 

Since, by hypotheses ^2) and .A3), vq G -D(^) and -jrt G D(y/L), the regularity of the limit 
function f in Theorem 2 can be precised further: it is the unique strong solution of Cauchy 
problem P-flTDl) and satisfies (|23|). 
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